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SLENDER  BODY  THEORY  FOR  THREE-DIMENSIONAL  BOUNDARY  LAYER 

t 

INDUCED-  POTENTIAL  FLOWS 

t  * 

By  Stanley  G.  Rubin  and  Frank  J.  Mummolo 

Polytechnic  Institute  of  Brooklyn,  Farmingdale,  New  York 

ABSTRACT 

The  application  of  slender  body  theory  for  the  evaluation  of 
three-dimensional  boundary  layer  induced  surface  crossflow  and 
streamwise  velocites  is  considered.  The  method  is  applicable  to 
sub-  and  supersonic  flows  when  small  perturbation  theory  appl3.es, 
and  when  the  Reynolds  number  is  large  so  that  the  thin  boundary 
layer  approximation  is  valid.  The  resulting  potential  problem  is, 
in  general,  reduced  to  a  two-dimensional  consideration  of  the  flow 
over  an  expanding  cylinder  with  porous  boundary  conditions.  As 
a  model  problem,  the  induced  surface  velocities  are  determined  for 
the  constant  density  flow  over  a  body  with  constant  elliptical 
cross-section.  The  limiting  solutions  for  a  flat  plate  of  finite 
span  and  a  near  circular  cross-section  are  obtained  in  a  simple 
analytic  form.  In  the  former  case,  within  the  limitations  of 
slender  body  theory,  the  results  are  in  exact  agr  ;.ement  with  the 
complete  three-dimensional  solution  for  this  geometry. 
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i.  introduction 

Three-dimensional  boundary  layer  induced  crossflow  analyses, 

in  particular  on  surfaces  with  locally  large  or  even  infinite 

* 

transverse  curvature,  belong  to  a  class  of  problems  for  which  the 
current  state-of-the-art  is  still  inadequate.  Only  a  limited  group 
of  analyses,  for  very  specific  geometries,  is  presently  available. 

These  include  the  quarter  plate  investigations  of  Stewartson  and 
Howarth  (1960)  and  Stewartson  (1961),  and  the  cruciform  studies 
of  Rubin  (1966),  Pal  and  Rubin  (1971)  and  Rubin  and  Grossman  (1971). 

One  of  the  difficulties  associated  with  the  determination  of 
boundary  layer  interactions  of  this  type  is  the  complexity  of  the 
three-dimensional  potential  problem  that  arises  when  the  displacement 
induced  flow  field  is  to  be  evaluated.  From  this  potential  flow 
solution,  boundary  layer  induced  streamwise  velocities  and  cross- 
flows  are  determined.  With  a  systematic  application  of  singular 
perturbation  methodology,  the  induced  velocity  distributions  along 
the  surface  become  the  necessary  asymptotic  matching  conditions 
for  the  crossflow  boundary  layer  analysis,  as  well  as  for  higher 
order  boundary  layer  considerations.  Special  care  must  be  exercised 
at  points  of  large  transverse  curvature  such  as  the  edge  or  corner 
intersection  line.  In  these  boundary  regions,  crossflow  diffusion 
is  important  even  for  the  lowest  order  theory. 

Since  the  crossflow  and  second-order  boundary  layer  analyses 
require,  as  asymptotic  boundary  conditions,  only  the  potential 
velocities  evaluated  along  the  surface,  the  slender  body  theory 
of  Munk,  Jones  and  Ward  (M-J-W) ,  for  many  years  a  powerful 

★ 

i.e.,  when  the  local  cross-sectional  curvature  multiplied  by  the  local 
boundary  layer  thickness  is  0(1)  or  larger. 


approximation  in  -'nviscid  aerodynamic  theory,  is  adapted  here  for 
application  to  high  Reynolds  number  viscous  flows. 

When  considering  parallel  flows  over  bodies  with  locally  small 
longitudinal  curvature,  the  concepts  of  small  perturbation  theory 
apply  (Sears  (1960) ) ,  The  governing  equations  in  terms  of  a 
velocity  potential  cp  reduce  to  the  simplified  linear  Prandtl- 
Glauert  equation  for  sub-  or  supersonic  flow; 

S2cp+cp+cp=0  (1) 

where  s  denotes  the  stream  direction,  y  and  z  the  cross  plane 
directions  (see,  figure  1);  R^l-M^;  M^  is  the  stream  Mach  number. 
For  incompressible  flow  p3=l  and  (1)  is  exact. 

For  slender  bodies  with  slowly  varying  cross-sectional  areas, 
the  M-J-W  theory  specifies  that  in  the  vicinity  of  the  surface,  a 
first  approximation  for  (1)  reduces  to 

cPyy  +  C£>2Z  =  °  (2) 

so  that  the  solution  of  (1)  is  of  the  form 

cp(s,y,z)  =  cp(y,z;s)  +bQ(s)  (3) 


where 


b  (S)  =  r  §4 

Vs;  4tt  in  4s { £-s)  4n  JQ  | 

b  (s).S^Slln  fS 

o  2tt  2s  2tt  J  s-x 


1  S/(x)-S/(s) 


ix  (M  <  1)  (4a) 


(M  >  1) 

CO 


S(s)  denotes  the  local  cross-sectional  area  of  a  body  of 
length  l,  see  Ward  (1955) . 

Therefore,  the  original  three-dimensional  potential  problem 
(1)  is  replaced  by  a  much  simpler  two-dimensional  problem  (2)  for 
the  induced  flow  near  the  surface.  This  remarkable  simplification 
allows  for  the  solution  of  a  wide  variety  of  aerodynamic  problems 
for  which  the  general  potential  solution  can  only  be  represented 
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by  quadratures. 

In  view  of  the  fact  that"  it  is  precisely  the  surface  velocities 
that  are  required  for  evaluating  the  induced  cross  flow  and  higher 
order  boundary  layers,  the  M-J-W  procedure  lends  itself  to  the 
evaluation  of  the  boundary  layer  induced  potential  flow.  A 
complete  discussion  of  this  procedure  as  well  as  application  to  an 
elliptic  cylinder  is  presented  in  this  paper.  Consideration  of 
specific  cross  flow  boundary  layers  as  well  as  the  general  effects 
of  large  transverse  curvature  in  boundary  layer  theory  will  be  the 
subject  of  future  considerations. 

2 .  Formulation 

Consider  the  uniform  flow  over  a  cylindrical  body  of  arbitrary 
cross  section  at  zero  incidence.  Local  body  slopes  are  assumed 
to  be  small  so  that  the  uniform  flow  is  only  slightly  disturbed 
by  the  presence  of  the  body  and  small  perturbation  theory  as 
prescribed  by  (1)  applies.  Furthermore,  variations  in  cross-sectional 
area  are  also  assumed  to  be  small  so  that  near  the  surfa.ce  the 
M-J-W  theory  with  (2)  and  (3)  apply  for  the  inviscid  potential 
problem. 

The  Reynolds  number  Rg=Ua/v  based  on  freestream  values  and 
a  typical  body  dimension  a  is  assumed  to  be  large  so  that  a  thin 
boundary  layer  develops  over  the  body.  Therefore,  the  boundary 
layer  induced  potential  flow  is  also  governed  by  (1)  and  in  view 
>f  the  linearity  of  (1)  and  the  small  disturbance  boundary  conditions, 

*The  effects  of  small  angles  of  attach  or  yaw,  while  m  f  precluded 
by  the  analysis,  are  not  specifically  considered  in  this  paper. 
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Sears  (1960) ,  the  viscous  and  geometric  problems  can  be  treated 

-k 

independently  to  first-order  in  Re  and  t,  a  characteristic 
thickness  ratio. 

The  geometric  inviscid  problem  has  been  discussed  in  great 
detail,  see  Ward  (1955),  so  that  we  shall  confine  our  remarks  to 
the  boundary  layer  induced  flow  over  cylinders  of  constant  cross- 
sectional  area,  t=0.  For  these  geometries,  provided  that  the  trans¬ 
verse  curvature  of  the  cylinder  is  nowhere  large,  a  Blasius  boundary 

r 

layer  with  thickness  6-s^2  forms  along  the  generators  of  the  cylinder, 
Cooke  (1957) .  Regions  with  locally  large  or  infinite  transverse 
curvature  (boundary  regions),  or  the  flow  near  a  blunted  nose  cap 
attached  to  the  cylinder,  must  be  treated  separately  by  local 
singular  perturbation  analyses,  see  Rubin  (1966)  and  Van  Dyke  (1964) . 
For  the  downstream  flow  considered  here,  the  shape  of  the  nose  is 
immaterial,  see  Seban  and  Bond  (1951) . 

It  should  be  emphasized  that  the  method  to  be  presented  herein 
requires  only  that  viscous  displacement  effects  are  small,  i.e.,  thin 
boundary  layers.  The  analysis  is  applicable  to  geometries  with  non¬ 
constant  cross-sections  provided  the  longitudinal  curvature  is  small, 
T«l.  For  these  cases  it  is  assumed  that  the  boundary  layer  solu¬ 
tion  is  determinable.  For  slender  bodies  the  inviscid  perturbation 
streamwise  velocity  is  of  order  T3lnT.  If  we  accept  an  error  of 


this  order  for  the  boundary  layer  formulation,  the  lowest-order 
solution  always  reduces  to  the  Blasius  function. 


The  surface  boundary  layer  perturbs  the  outer  potential  flow 
in  much  the  same  way  as  streamwise  variations  in  body  cross-sectional 
area.  The  viscous  displacement  induces  crossflows  and  higher  order 
streamwise  velocities.  When  evaluated  at  the  surface,  these 
velocities  serve  as  asymptotic  matching  conditions  for  the  boundory 
layer  determinations.  The  effective  displacement  body  remains 
slender,  so  that  the  M-J-W  theory  with  (2-4)  should  apply  near 
the  surface.  S(s)  can  now  be  interpreted  as  the  modified  body 
cross-section,  taken  to  include  the  displacement  boundary  layer. 

For  a  general  slender  body  at  zero  incidence  the  familiar 
inviscid  M-J-W  slender  body  formulation  reduces  to  a  two-dimensional 
problem  (2)  for  the  flow  over  an  expanding  body.  The  boundary 
condition  applied  at  the  surface  relates  the  normal  gradient  of  the 
potential  function  cp  to  the  local  streamwise  body  slope.  For  the 
boundary  layer  induced  potential  flow,  the  slender  body  theory  also 
leads  to  a  two-dimensional  problem  (2)  for  the  flow  over  an 
expanding  body.  The  boundary  values  are  obtained  by  asymptotic 
matching  to  the  inner  boundary  layer.  For  large  Re  this  condition 
is  applied  at  the  surface  of  the  cylinder  (see  footnote  on  page  8) . 

In  view  of  vhe  linearity  of  this  potential  problem,  the  complete 
solution  can  be  obtained  by  superposition  of  the  geometric  and 
viscous  induced  flows. 

The  details  of  the  analysis  aie  presented  in  the  next  section 
where  the  incompressible  viscous  flow  over  a  cylinder  v>'ith  constant 
elliptic  cross-section  is  considered.  The  complete  solution  is 


6 


described  and  the  axisymmetric  and  zero  thickness  limits  are 
obtained.  A  particularly  simple  form  of  the  solution  results  for 
cross-sections  that  are  nearly  circular.  The  zero  thickness  limit 
is  compared  with  the  three-dimensional  solution  obtained  from 

an  extension  of  the  quarter-plate  analysis  of  Stewartson  and 
Howarth  (1960) . 

3.  Solution  For  An  Elliptic  Cylinder 

Consider  the  incompressible  flow  over  a  cylinder  of  constant 

elliptic  cross~section .  The  generators  of  the  cylinder  are  taken 

parallel  to  the  freestream  with  l >  s >0  (figure  1).  The  Reynolds 

number  based  on  the  semi-major  axis  is  large  so  that  boundary  layer 

theory  applies  in  the  usual  manner.  It  is  assumed  here  that  the 

transverse  curvature  of  the  cr^ss-section  is  small  with  respect 

to  the  boundary  layer  thic-Kness,  although  in  a  more  general 

★ 

consideration  this  audition  may  be  violated  locally  along  discrete 
generators.  This  would  occur  at  the  ends  of  a  highly  eccentric 
elliptic  cross-section  or  in  a  l_miting  case  at  the  edges  of  a 
plate  of  zero  thickness.  These  singular  lines  diffuse  into  what 
are  termed  boundary  regions  and  must  be  treated  separately  (see 
Stewartson  (1961)  or  Rubin  (1966) )  .  To  lowest-order,  however „  the 
viscous  boundary  regions  do  not  affect  the  flow  behavior  in  the 
adjoining  three-dimensional  boundary  layers  or  the  outer  inviscid 
potential  flow. 

As  stated  previously,  to  lowest-order  the  boundary  layer  on 
a  constant  area  cylinder,  downstream  of  a  leading  edge  or  nose 

*  In  small  local  regions  that  are  of  zero  thickness  as  v-0 . 


'  ~‘,7TTZ?iW^\’p?Ar 


region  and  away  from  boundary  regions,  is  identical  with  the 


constant  pressure  flat  plate  boundary  layer.  Transverse  curvature 


enters  only  in  higher  order  approximations.  This  solution,  due 


to  Biasius,  is  well-known,  Van  Dyke  (1964) .  It  is  sufficient  to 


note  here  that  the  viscous  induced  displacement  velocity  v  directed 


normal  to  the  surface  has  the  following  asymptotic  behavior  for 


large  values  of  the  boundary  layer  variable  v^=y(u/2vs)  2: 


lim  v(  s,  v) 

'n— *03  U 


0.860 


and  the  displacement  thickness  A  is  given  ?jy  the  relation 


A*  =  1.72s/Rs^ 


where  Rs=Us/v?  U  denotes  the  undisturbed  stream  and  (y,x, s)  are 


the  surface  normal,  azimuthal  and  streamwise  coordinates, 


respectively,  defined  in  figure  1. 


For  large  Reynolds  numbers,  R  »1,  the  inviscid  outer  flow 

s 


is  only  slightly  perturbed  by  the  formation  of  the  thin  boundary 


layer  on  the  surface  and  therefore  the  solution  for  the  incompressible 


displacement  induced  potential  flow  can  be  determined  from  ('.) 


with  M  =0 . 


A  far  field  velocity  decay  condition 


7  co-0 


applies  at  large  distances  from  the  body  and  a  smooth  match  with 


the  inner  boundary  layer  must  be  prescribed.  This  surface  matching 


condition  is  specified  by 


lim  n  •  \’cc( s,y,  x)  =.lim  v/U  ~  0.860^  ^ 

A  .  _  S 


When  the  largest  cross-section  dimension  is  of  die  cader  of  the  displacement 
thickness,  the  boundary  condition  cannot  be  app? ied  at  y=0,  but  only 
at  y=A  •  This  point  will  be  illustrated  in  a  later  examole. 


t »J^9Wi^}Wn® ' jvV is>3 <r ^>*’~  -.--**  ,  • 


In  most  cases  the  boundary  condition  (7b)  can  be  imposed 


directly  on  the  surface  y=0;  .however,  the  resulting  three-dimensional 


potential  problem  is  still  a  formidable  one.  This  is  particularly 


true  with  complex  body  cross-sections.  To  the  authors'  knowledge, 


three-dimensional  boundary  layer  induced  potential  flows  have  been 


calculated  only  in  a  few  cases  with  simple  geometries,  e  g.,  the 


quarter-plate  of  Stewartson  (1961)  and  the  cruciform 


of  Rubin  (1966) .  In  the  former  case  only  integrated  surface 


conditions  were  obtained. 


In  order  to  generate  surface  values  for  application  in  higher 


order  boundary  iayei  theory  and  to  evaluate  the  induced  surface 


cross  flow,  the  three-dimensional  potential  problem  can  be  reduced 


to  a  qua si- two-dimensional  form  with  the  M-J-W  slender  body 
approximation.  The  displacement  body  prescribed  by  (6)  describes  a 


slender  body,  and  the  governing  equation  (2)  with  (3)  and  (4) 


applies.  As  stated  previously,  S(s)  denotes  the  cross-sectional 


area  of  the  body  plus  displacement  thickness. 


In  general,  the  solution  of  (2)  satisfying  (7)  can  be 


investigated  with  conformal  mapping  techniques.  However,  for  the 


elliptic  :ross~section  considered  here,  it  is  convenient  to  trans¬ 


form  to  an  elliptic  coordinate  frame  in  which  tne  general  solution 


is  of  the  form 


cp  —  a  (s)  f  +  c  ( s)  +  2  a  (s)e  ^cosnf) 
°  n=l  n 


if  „ 

It.  should  be  noted  that  with  S'(s)~sp,  the  integral  in  (4a)  is 
divergent  for  semi-infinite  bodies  (t-®)  with  p>0;  however,  the 
perturbation  velocities  depend  on  b  '(s)  and  are  well-defined  for 

p>0 .  ° 

9 
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where  (§,9)  are  elliptic  coordinates,  figure  2,  defined  by 

z=ccosh  ^cos'i  (9a) 

y=csinh|3in0  (9b) 

r2=y2  +  z2  ;  ju=tan"  L  (y/z)  (9c) 

§  =  In  (— )  +  ^ln  jl+  [  (1-  p-)-'J  +4  (—)  3  sin2  d]^ 

+  rS  1-  (pr)cos2'J,  +  {(1-  ~)  s+4  (p-)  sin2  i}'5)  3^  f  (9d) 

and  c3=a2-b2,  e=c/a= (cosh?^)  ,  where  a  and  b  are  the  semi-major 

and  minor  axes,  respectively,  and  e  is  the  eccentricity. 

The  boundary  condition  (7b)  becomes 

lim  [c  (cosh3 --cos3  9)  ^]-1  ~ ,  9;  s)  =  (0.860/R  ^)  (a/s)^= 

-  7  "  o  5  ‘  ^  as 

J°  (10) 

where  the  surface  is  denoted  by  the  ellipse  ?=*Q. 

The  coefficient  a  are  defined  by 

n 


>0-£f  ’'ff(V?!s)d0 


V  I  ff<?0.=:s)cosn,d8 

o 


(11a) 


(11b) 


where  -^f(?  ,9,s)  is  determined  from  (10). 
o '  "O 

The  coefficient  2ra  can  be  shown  to  represent  the  derivative  of 

o 

the  "effective"  cross-sectional  surface  area  of  the  ellipse  denoted 
here  as  S'(s).  This  result  is  exact  when  the  boundary  condition 
(7b)  is  satisfied  at  the  displacement  surface  and  in  certain 
limiting  cases  where  the  boundary  condition  is  transferred  to  the 
surface  (10),  e.g.,  a  circular  cross  section  or  flat  plate.  In 
general,  when  the  boundary  condition  is  transferred  to  the  surface. 


This  result  can  easily  be  proven  by  a  mass  conservation  argument. 


"e&sasBEP sg&gm.  ^3^  3«l  - 


as  given  by  (10) ,  a^  differs  from  S'(s;  by  a  term  of  order  A  2. 
This  difference  is  negligible  in  the  present  analysis  but  would 
reappear  if  higher  order  solutions  for  the  outer  inviscid  flow 
were  considered.  Therefore,  the  surface  area  terms  S'(s)  in  (4) 


are  represented  by 


S ' (s)  =2rrao=  { lim  J  c  (cosh2  ?-cos2  0)^30}--^- 


F-S  O 
•O 


S/(s)=4a  E (~, e) 

2  as 


(12a) 


where  E  is  a  complete  elliptic  integral  of  the  second  kind,  defined 
in  Gradshteyn  and  Ryzhik  (1965) ,  p.  909. 

Some  limiting  cases  are 

(1)  For  a  circle  of  cadius  a:  c-0,  ccosh^  -a,  so  that 

*  ® 


S/(s)=2TTa 

ds 

* 

(2)  For  a  circle  of  radius  a=0(A  ):  c-0, 

* 

lim  ccosh^-a+A  and  (12)  gives 

5-S0  *  3  A* 

°  S'(s)=2Tr(a+A  )  f” 


(12b) 


(12c) 


This  case  is  not  physically  meaningful,  as  the  ratio  a/A  =0(1) 
ano  the  Blasius  solution  is  no  longer  appropriate;  however,  it 
provides  an  interesting  limit  when  a-0  because  the  body  is  then 
reduced  to  a  paraboloid  of  revolution.  The  exact  solution  for  this 

problem  was  aiven  bv  Van  Dvke  (1964^ . 


(3)  For  a  flat  plate  of  span  2a:  b-0,  c-a,  Z  -0  and 

o 


(12a)  gives 


S'(s)  =4a 


(12d) 


The  coefficient  cq(s)  in  (8)  is  defined  so  that  the  expressions 
(4)  obtained  bv  Ward  (1955)  are  applicable  in  that  form,  see 


Kahane  and  Solarski  (1953) . 


co(s)  =  ao(s)  In  - 


(12e) 


Therefore,  with  (4),  (10-12e),and  notinq  that  odd  coefficients  a 

r 

vanish,  the  solution  for  o  becomes 

c?(s,§,G)  =  — ^  U+ln  -f)+b  (s)  + 

2iT  2  O 


-1.72c  00  e  "o'  \  r‘v  *■  P  > 

~~T~  £  ~  (  (cosh  f  —cos®  6)  2cos  2n  £  dQ  ( cos  2n  6  (13) 

rrR  ^  n=l  n  <  o  °  ! 


“2n(?-0  Un/2 


The  induced  velocities  at  the  surface  of  the  ellipse  ?=§  are 

f  ■  V 

a 


_£  — -~n 9  { v/l-e®cosa3  cos2n-dS(] 
n=l  o 


(14a) 


IT  =  3  U-e*cos3  6)  _liCg  = 


3. 44. a,  2  ®  sin2nS  0  ,  ,  h 

- 2  ~ — - — r  (l-e3cos3S  2con2n3dp  (14b) 

nRg  2  s  r,=i  ( 1-e®  cos®  0)  2  vq 

Expressions  (14a)  and  (14b)  represent  the  boundary  layer  induced 
velocities  evaluated  at  the  surface  of  a  body  having  a  constant 
elliptic  cross-section.  The  crossflow  boundary  layer  and  first 
order  correction  to  the  streamwise  Blasius  profile  require  these 
values  as  outer  matching  conditions.  Previous  boundary  layer 
analyses  of  transverse  curvature  effects  on  cylinders  of  circular 
and  more  general  cross-section  have  erroneously  neglected  these 
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Cooke  (1957), 


displacement  induced  boundary  conditions,  e.g., 

Seban  and  Bond  (1951).  a  dis'cussion  of  the  errors  arising  from 
the  neglect  of  these  displacement  effects  for  the  case  of  a  circular 
cylinder  was  given  by  Van  Dyke  (197*).  The  complete  boundary 
layer  analysis  is  currently  under  investigation.  For  the  remainder 
of  the  present  paper,  the  solutions  (14)  for  the  induced  velocities 
will  be  discussed  for  different  eccentricities;  the  limiting  circular 
and  flat  plate  results  will  be  obtained  for  comparison  with  known 
exact  solutions.  This  will  establish  the  validity  of  the  M-J-W 
theory  to  boundary  layer  problems  and  serve  as  a  model  for  examinations 
of  more  complex  cross-sections  where  exact  three-dimensior.ol  solu¬ 
tions  are  not  presently  available. 

4.  Flat  Plate  of  Finite  Span 

The  flat  plate  represents  the  simplest  geometry  for  which 
displacement  induced  crossflows  appear.  With  the  formulation 
presented  herein,  the  surface  velocities  can  be  evaluated  from  the 

general  solution  of  section  3  by  a  direct  limiting  process  ?  -0, 

o 

c-a.  The  solution  for  this  geometry  is  also  easily  obtainable  by 
the  appropriate  application  of  Green's  functions j  two-dimensional 
for  the  slender  body  theory  discussed  here  and  three-dim'’  is*onal  for 
the  classical  formulation.  This  latter  approach  has  been  considered 
by  Stewartson  and  Howartn  (1960)  for  a  quarter  -infinite  flat  plate 
and  is  readily  extended  to  the  finite-span  case. 
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A.  Three-Dimensional  Solution 

Consider  a  plate  of  span  2a  as  shown  in  figure  3.  Away 
from  the  side  edg^  r^undary  regions,  a  Blasius  boundary  layer  forms 
on  the  plate  and  the  three-dimensional  displacement  induced 
potential  problem  becomes 

CD  +  CO  +  CD  =0 

vs  s  y  y  z  z 


with 


vcp-o  as  (y?+22)  2  -  oo 


-U  .  . 

and  (cPy)y=0=0.860URs  sgny;  1 2  | <a ,  s>0 

and  (cpy)y=0  =  0  elsewhere- 

Following  Stewartson  and  Howarth  (1960),  the  solution  is  given 


,  x  0.860U  r  dx  r 

-2rrcp(s,y,z)  =  - r~  j  ~T  j 

D  r.  ^  -  V  2 


R„  2  o  X2  -a  [  ( r-z)  "+y2+  (r-s) 


It  follows  that 

oo  a 

’  _  -0.860  iv  j  dX  P 

▼  -  O  Al  r  r  "l  :  ' 


860  J  V  I  dX  |'  _ (  §-z)d  c _ 

11  U  l'o  ‘'-a  r(?-z) -+y2+(X-s) 3 


_ oo  a 

u  -0.860  iv  (X-s)dX 


'o  X 


-a  r  (  ~,-z)  ?+y3+  (X-s)  c 


On  the  plate  y=0,  with  a/s<"<l,  the  crossflow  velocity  takes  the  form 


,w,  -0.860  r a  —  z i 

'u’y-0  =  109  fe- 


(15a) 


Along  the  plate  centerline,  y=0,  z=0  and  with  a/s«l,  the  streamwise 
perturbation  becomes 


(15b) 


,u4  _  -0.860  ra  y,  a  .  ,  »  , 

U  y=0=z  -  -rr?(l0«  8?  +  3 

s 


while  at  the  plate  edge  we  obtain 


,.ux  _  -0,860,3^3/2,..  ,  a  > 

u  y-o  -  [i09(4S>  t 

_  „  TTK 

z=a  a 


(15c) 


B.  Flat  Plate  as  the  Limit  of  the  Elliptic  Cylinder  with 
M-J-W  Theory 

This  limit  can  be  obtained  from  (14)  by  setting  5qs0 
and  noting  that  from  (12d)  that 

a  fs)  -  islL.  S/.(s) 

V®'  "  TT  CIS  "  2TT 

The  streamwije  perturbation  velocity  (14a) ,  on  the  plate  centerline 
9=rr/2  becomes 

.u.  S^r.  a  2.  (-1)  ->  ,  .  ,  .  * 

Vy=0=z_  2rr  ln  2  +S.  n(4na-l)-1  +  bo  (s) 

n=l 

= -^-5--  <-)  Clog  ^  +  1]  (16) 

s 

and  at  the  edge  (15c)  is  recovered.  See  Gradshteyn  and  Ryzhik  (1965) , 


p.  46. 


The  crossflow  on  the  plate  is  given  by 

(—)  =  ~s ' (s)  2  sin2n9 

'U  y=0  ”  TTa|sin9|  ,  4na-l 

n=l 


Here  b^(s)  denotes  the  derivative  of  (4a)  in  the  iimit  t-co  . 


It  is  easily  shown  that 


4 

I  sin0 


®  sin2n  8 
S.  4n8  -1 


l°g  (77?) 
a+z 


where  from  ^9b)  z=acose.  Therefore 


,w,  _  -0.860 

%>y=0  -  r 


l09  <|ff) 


These  results  (16)  ana  (17)  are  in  exact  agreement  with 
three-dimensional  solutions  (15)  . 

C.  M'-J-W  Theory  for  the  Flat  Plate.  Direct  Formulation 


The  slender  body  solution  for  the  induced  velocities  on 
a  flat  plate  can  also  be  determined  by  direct  application  of  a 
two-dimensional  Green's  function.  With  the  governing  equation 


for  cp(y ,  z  ;  s ) 


A  ,  A  _  _ 
CO  +  CD  ~  0 

Yyy  ^zz 


and  the  boundary  conditions 


A 

vcp  -  0 


as  (y s+z2 )  -  oo 


, a  i  0 .860  i  i 

(c^)y=0  =  — sgn  y  |z|  < 
Rs 

(C^)y=0  *•  0  |z|  >  a 


it  follows  that 


Therefore, 


o_a  0.860  i*  .  r  a  , 

2  TTcp  =  - r—  :  log[  (X-z)  e+y2  ]dX 

Rs1 


w  -0.860  „ 

u =  mr  1091 


ISZillixL 

(a+zr+r 


’^.i’i'.  ,!H*S3S!a:s=Trt.’fR' 


and 

{^z=0=  (f}  tlog ( aa+y 3 ) ^-1-  J  tan"1  J]  +  b'(s)  (19) 

ttr  y  u 

s 

On  the  surface  y=0  the  solutions  (15)  are  recovered  for  u/U  and 
w/U  respectively. 

5.  Quasi-Axisymmetric  Flow 

For  geometries  where  the  cross-sections  are  close  to  circular 
(e«l)  ,  the  leading  terms  in  an  expansion  of  (14)  for  small  e  take 
a  particularly  simple  form.  W c  find  that 
-0.430  ,a, 3/2, 


and 


u. 

Ra 


f  (P  [  (1"‘ 4  - »  PJ ~ >  ln  +2  +“COs2  9]+0  (e4  )  (18a) 


O^^^^sin29  +  °<e4> 


(18b) 


where  from  (12) 


S ' (s)  =  ■x--^"-a(l-~)  +  0(e4) 
R  *2  4 


The  relation 

tt/2 


I 


2m  „  , 

cos  x  cos2nx  dx  = 


0 

TT 


2m+l  vm-n 


for  m<n 

,2m  .  _ 

(  )  for  m>n 


has  been  required  in  obtaining  (18)  and  is  useful  if  additional 

terms  in  the  expansion  are  desired.  The  series  in  (14)  reduce  to 

2  n 

n  terms  for  the  e  coefficient.  The  expressions  (18)  are  presented 
here  as  limiting  velocities,  as  e-0,  for  the  general  solutions  with 
elliptic  cross-section  to  be  discussed  in  the  following  section. 

The  axisymmetric  solutions  (e=0)  have  .previously  been  obtained  by 
Van  Dyke  (1970) ,  who  also  discussed  the  second  order  boundary  layer 
and  corrected  an  earlier  paper  by  seban  and  Bond  (1951). 
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Finally,  the  special  limiting  case  of  a  needle  a-0,  e-0  leads 
to  the  flow  over  a  slender  paraboloid  of  revolution.  While  the 
Blasius  boundary  layer  is  incorrect  in  this  limit,  (see  Glauert 
and  Lighthill  (1955)),  the  exact  solution  for  a  paraboloid  is 
Known  and  provides  an  additional  check  on  the  theory.  Van  Dyke  (1964) . 
With  the  cylindrical  radius  defined  by 

Rq(s)  =  (2  e3s ) ^ 


Van  Dyke  has  shown  that 


cp  =  -y  log(r2/2e8s) 


=  -eV2s  (19) 

From  (4)  and  (14a),  with  (12b),  we  obtain 

jj  =  b^(s)  =  -  (0.660)  3Rg"*1  (20) 

With  the  cylindrical  radius  set  equal  to  the  effective  body  radius 
or  displacement  thickness  A*  from  (6),  we  find  that  (19)  and  (20) 
are  identical. 


6.  Elliptic  Cylinder 

The  flow  properties  at  the  surface  of  a  cylinder  of  elliptic 
cross-section  with  arbitrary  eccentricity  e  have  been  determined 
by  numerical  evaluation  of  an  appropriate  number  of  terms  in  the 
series  solutions  (14).  The  integrals  in  (14)  can  be  expressed  by 
means  of  recursion  relations  in  terms  of  elliptic  integrals  of  the 
first  and  second  kind,  see  Gradshteyn  and  Ryzhik  (1965).  it  was 
more  convenient  to  evaluate  the  integrals  directly  by  numerical  inte¬ 
gration.  A  three -point  Simpson's  rule  having  a  nominal  error  of  order  h‘* 


was  used ; h  represents  the  grid  spacing.  As  n  increases,  the  integrand 


exhibits  rapid  oscillations  and  an  increased  number  of  mesh  points 
are  required,  in  order  to  maintain  a  fixed  error,  twelve  intervals 
were  always  located  between  the  zeros  of  cos2n0  so  that  24n+1  points 
result  for  n=N. 

The  number  of  terms  N  required  to  achieve  a  specified  degree 
of  accuracy  increases  sharply  as  e-1.  The  summation  was  terminated 
when  percentage  differences  in  the  sum  S  and  S+5  terms  were  less 

-iO 

than  10  .  The  singular  behavior  in  the  flat  plate  limit  and 

in  particular  near  the  side  edge  is  not  unexpected.  For  e=l  a 
logarithmic  singularity  in  the  crossflow  velocity  (14b)  appears 
at  the  edge,  while  for  e<l,  w-0.  Also,  the  asymptotic  form  of  the 
integrals  in  (14)  shows  that  e=l  is  a  singular  limit,  since 

llm  ['"vfcpwi  cos2n0de~HA(^’M  f°r  an>’  £inite  M>0  (e<1) 

n-°°  o  (4n  -l)"1  for  e=l 

The  asymptotic  form  for  e=l  is  also  obvious  from  the  exact 
solution  (15a) . 

The  numerical  solutions  for  the  crossflow  velocity  are 
depicted  in  figure  4  and  confirm  the  asymptotic  predictions.  Note 
the  large  increase  in  the  maximum  number  of  terms  N  required  for 
the  series  solution  as  e-1.  For  e=0.999,  the  solution  closely 
follows  the  flat  plate  logarithmic  behavior  before  falling  rapidly 
near  the  edge.  The  solution  exhibits  boundary  layer  behavior  near 
the  edge,  with  the  peak  crossflow  velocity  approaching  the  edge 
as  e-1.  The  maximum  number  of  terms  required  for  the  series 
summation  generally  occurs  at  a  location  corresponding  to  the  peak 


w  value.  A  boundary  region  forms  near  the  edge  and  the  three- 
dimensional  viscous  flow  must-be  reevaluated. 

Similar  behavior  for  the  streamwise  velocity  perturbation 
is  shown  in  figure  5.  While  the  streamwise  velocity  distribution 
is  continuous  as  e~l,  the  velocity  gradient  also  exhibits  a 
logarithmic  singularity.  Also  depicted  on  the  figures  4  and  5  are 
solutions  for  several  values  of  e  between  zero  and  one.  The 
simple  quasi-axisymmetric  expansions  for  e«l  are  given  by  (18a) 
and  (18b)  „-nd  are  shown  i-n  figures  C  and  7.  The  one  term  expansions 
agree  quite  well  with  the  numerical  solutions  for  e<0.25. 

7 .  Summary 

When  one  attempts  to  analyze  the  perturbation  velocities,  and 
in  particular  the  cross  flows,  induced  by  thin  boundary  layers  on 
three-dimensional  geometries,  a  difficult  potential  problem  results. 

To  date,  there  are  only  a  minimum  number  of  such  exact  solutions 
available.  When  considering  the  flow  over  thin  bodies,  it  1  as 
been  shown  herein  that  slender  body  theory,  for  many  years  a 
powerful  tc o 1  of  classical  aerodynamics,  can  be  adapted  for  the 
determination  of  boundary  layer  induced  flows.  In  essence,  the 
potential  problem  is  reduced  to  the  consideration  of  a  two-dimensional 
flow  over  an  expanding  body  with  porous  boundary  conditions.  This 
results  in  a  considerable  simplification. 

The  solutions  are  valid  near  the  body  surface  and  provide 
insight  into  the  types  of  crossflows  generated  by  three-dimensional 
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configurations,  as  well  as  supplying  asymptotic  matching  conditions 
for  higher-order  boundary  layer  theory.  These  displacement 
induced  velocity  boundary  conditions  have  generally  been  neglected 
in  previous  perturbation  analyses  dealing  with  transverse  curvature 
effects  on  boundary  layer  behavior. 

The  constant  density  flow  over  a  cylinder  of 
constant  elliptical  cross-section  has  been  considered.  Complete 
solutions  have  been  obtained  in  the  form  of  a  Fourier  series. 

For  the  limiting  cases  of  a  finite  span  flat  plate  or  near  circular 
geometry,  simplified  analytic  forms  have  been  obtained.  The 
surface  solutions  for  the  plate  are  in  perfect  agreement  with  the 
exact  three-dimensional  result  obtained  by  a  simple  extension  of 
the  method  of  Stewcrtson  and  Howarth  (1960) . 

In  the  general  case,  a  peak  cross  flow  velocity  is  predicted. 
The  location  of  this  maximum  moves  toward  the  side  region  of  the 
ellipse  as  the  eccentr icity  e->l.  In  this  limit,  the  crossflow 
exhibits  a  logarithmic  behavior,  corresponding  to  the  flat  plate 
solutions,  up  to  the  peak  location  and  then  falls  rapidly  to  zero 
at  the  symmetry  line.  For  e-1,  the  application  of  a  Blasius 
boundary  condition  near  the  edge  is  incorrect/  as  transverse 
curvature  is  quite  large  there.  Further  study  will  have  to  consider 
the  side  edge  flow  more  exactly  in  order  to  determine  the  precise 
nature  of  the  developing  boundary  region  where  curvature  and  cross¬ 
diffusion  effects  become  important. 
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FIG.  6  STREAMWSSE  PERTURBATION 


a/s  «  0.1 


FIG.  7  CROSS  FLOW:  Q  AS  THEORY 


